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A two-level quantum system coherently driven 
by a resonant electromagnetic field oscillates si- 
nusoidally between the two levels at frequency 
f2 which is proportional to the field amplitude 
This phenomenon, known as the Rabi oscilla- 
tion, has been at the heart of atomic, molecular 
and optical physics since the seminal work of its 
namesake and coauthors [2|. Notably, Rabi os- 
cillations in isolated single atoms or dilute gases 
form the basis for metrological applications such 
as atomic clocks and precision measurements of 
physical constants Both inhomogeneous dis- 
tribution of coupling strength to the field and in- 
teractions between individual atoms reduce the 
visibility of the oscillation and may even suppress 
it completely. A remarkable transformation takes 
place in the limit where only a single excitation 
can be present in the sample due to either initial 
conditions or atomic interactions: there arises a 
collective, many-body Rabi oscillation at a fre- 
quency y/~NQ involving all N ^> 1 atoms in the 
sample [4\. This is true even for inhomogeneous 
atom-field coupling distributions, where single- 
atom Rabi oscillations may be invisible. When 
one of the two levels is a strongly interacting Ry- 
dberg level, many-body Rabi oscillations emerge 
as a consequence of the Rydberg excitation block- 
ade. Lukin and coauthors outlined an approach to 
quantum information processing based on this ef- 
fect [5|. Here we report initial observations of co- 
herent many-body Rabi oscillations between the 
ground level and a Rydberg level using several 
hundred cold rubidium atoms. The strongly pro- 
nounced oscillations indicate a nearly complete 
excitation blockade of the entire mesoscopic en- 
semble by a single excited atom. The results pave 
the way towards quantum computation and sim- 
ulation using ensembles of atoms. 

A two-level quantum system coherently driven by a 
quasi-resonant electromagnetic field is one of the center- 
pieces of modern quantum physics (T), Q. A wide ar- 
ray of two-level systems have been realized, with atoms, 
molecules, nuclei, and Josephson junctions being some of 
the prominent settings. More than half a century ago 
Dicke recognized that an atomic ensemble coupled to an 
electromagnetic field cannot always be treated as a collec- 
tion of independent atoms [|[ . His ground-breaking work 
gave rise to a rich field of collective atom-field interaction 
physics [f|. 

A key prediction of Dicke's theory is that under cer- 



tain conditions atom-field coupling is enhanced by a fac- 
tor ~ \/N when compared to one atom. Collectively- 
enhanced atom-field coupling has since been observed in 
a variety of settings involving either the emission or ab- 
sorption of radiation. A coherent multi-atom Rabi os- 
cillation at a frequency y/~NQ is a particularly dramatic 
manifestation of quantum mechanics at work on meso- 
scopic scales, where an entire ensemble exhibits the dy- 
namical behavior of a single two-level system. In 2001, 
Lukin et al. proposed to realize many-body Rabi oscil- 
lations in ensembles of atoms driven by a laser tuned to 
a Rydberg level, and outlined designs for scalable quan- 
tum gates for quantum computation and simulation and 
generation of entangled collective states for metrology 
beyond the standard quantum limit. 

When an atom is promoted into a Rydberg level with 
principal quantum number n, the valence electron is in 
an orbit that is ~ n 2 larger than that of the ground- 
level atom. The atomic dipole moment is correspond- 
ingly larger, so that the dipole-dipole interaction of two 
atoms is increased by ~ n 4 Q. For n ~ 100 the inter- 
actions are sufficiently strong that for two atoms sepa- 
rated by a distance ~ 10 fim the associated energy shift 
may prevent the second atom from being excited. The 
excitation blockade gives rise to an oscillation between 
the collective ground state \G) = YiiLi an d the state 
\R) = l/v^Eili lff)i-k>i-lfl)jv in which one of the N 
atoms is in the Rydberg level |r), with frequency y/~NQ 
[E S @-Q3 • The average number (N) r of atoms in level 
\r) is given by: 
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(l/VEr=i fi f)Er=i^l5)i-|r)i...|p)iV ! see Meth- 
ods section. For two atoms the Rydberg blockade 
and the accompanying v2 enhancement of the Rabi 
oscillation frequency have been observed by Urban et al. 
11 and Gaetan et al. [l2|. Over the past decade there 
have been numerous studies of the many-atom Rydberg 
blockade, however none of them, without exception, 
demonstrated either a blockade by a single atom or the 
many-body Rabi flopping (lljl - Eol ] . 

Here we report observations of many-body Rabi oscil- 
lations of a mesoscopic (a ~ 15 fim) ensemble of rubidium 
atoms in the regime of the Rydberg excitation blockade 
by just one atom. In order to achieve this, the interaction 
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FIG. 1: Rydberg excitation of an atomic ensemble, a, Illustration of the excitation blockade of more than one Rydberg atom 
in the ensemble. The coherent laser driving couples the collective ground state |G) and the state with one Rydberg atom \R) with 
Rabi frequency -i/Nil. The doubly excited states \RR) are shifted in energy out of laser resonance by the strong atomic interactions, 
b, Single-atom energy levels for 87 Rb. Electronic, hyperfine, and Zeeman quantum numbers are shown. The detuning from the 
intermediate |5pi/ 2 ) level is Ai = —40 MHz. c, Probability P of photoelectric detection event per trial as a function of two-photon 
detuning A2 for level |102s 1 / 2 ). It shows the two rrij = ±1/2 Zeeman components split by the bias magnetic field. The solid curve 
is a sum of two Lorentzian functions fit with the 0.9 MHz widths (fwhm) of the peaks determined by the 1 /is excitation duration. 
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FIG. 2: Coherent many-body Rabi oscillations of a mesoscopic atomic ensemble. In panels (a)-(d), probability of photoelectric 
detection P as a function of the single-atom Rabi angle 6 is shown; upper level is |l02s!/ 2 ), excitation duration is r = 1 fxs. 

Solid curves are fits of the form P — ^Ae~ aB (1 — e~ be cos(y/ N e 0)) , see Methods section. The fit parameters (A, a, b, N e ) are: 
(4.3, 1.43, 1.70, 456) for a, (4.44, 1.43, 1.77, 397) for b, (3.24, 1.14, 0.72, 243) for c and (2.56, 0.79, 0.86, 148) for d. (e) yfNl as a 
function of number of atoms N a determined from fluorescence measurements. The data are fit with a function Cy/N a , with the 
best-fit value C = 0.74. The inset shows a collective Bloch vector tipped by the angle \/N e 9 on the unit sphere corresponding to 
the many-atom states \G) and \R). The error bars represent ± one standard deviation (vM) for M photoelectric counting events. 



strength AEb = A^- (a) between a pair of atoms at a dis- 
tance equal to the ensemble size a must be greater than 
the spectral width 5u> of the exciting laser held. The du- 
ration of coherent atom-light interaction is limited to < 2 
/is by the finite coherence time of the ground-Rydberg 
transition caused by atomic motion [2fj| . Therefore, we 



choose a high-lying Rydberg level |r) = IIO2.S1/2), for 
which AEb ^ 5 MHz is sufficiently large. 

A gas of 87 Rb atoms of temperature T ~ 10 fiK and 
of peak density po — 10 12 cm -3 is prepared in an op- 
tical lattice. The lattice is shut off, and the atoms are 
driven in resonance between the ground \g) = |5si/2) 
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FIG. 3: Probability of photoelectric detection P and second- 
order intensity correlation function at zero delay g (2) (0) as 
a function of the single-atom Rabi angle 6. Excitation du- 
ration is t = 1 fxs and upper level is 1 102s 1( / 2 ). In panel a 
the solid curve is a fit as in Figure 2(a-d). The fit parameters 
(A,a,b,N e ) are (3.80, 1.48, 1.86,492). The error bars represent 
± one standard deviation (VM) for M photoelectric counting 
events. 
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FIG. 4: Probability of photoelectric detection P as a func- 
tion of the single-atom Rabi angle 6. Excitation duration is 
t = 1 /is. The solid curves are fits as in Figure 2(a-d), where the 
fit parameters (A,a,b,N e ) are (4.10,2.00,3.52,441) for n=90 
in a and (3.42,1.62,6.70,335) for n=81 in b, respectively. The 
error bars represent ± one standard deviation (\/M) for M pho- 
toelectric counting events. 



and a Rydberg |r) level with the two-photon Rabi fre- 
quency f2(r) = fii(r)02(r)/(2Ai) for a duration r = 1 
/is, with the corresponding single-atom excitation pulse 
area 9 = fi(0)r, Figure 1 (a) and (b). The transverse size 
(Gaussian waists w x ~ w y ~ 6 p,m) of the Rydberg exci- 
tation region is determined by the overlap of the nearly 
counter-propagating two-photon excitation laser fields Oi 
at 795 nm and ^2 at 475 nm. The longitudinal extent of 
the ensemble is determined by the sample size of waist 
w z w 11 /im along z. 

The measurement of the population of state \r) is 
achieved by the quantum state transfer onto a retrieved 
light held using a 1 /is long read-out field f?3 at 475 
nm, in resonance with the |102s!/ 2 ) ^ |5pi/2) transition 
[2TI [22j . The retrieved field is coupled into a single mode 
fiber followed by a beam splitter and a pair of single- 
photon detectors Di and D2. Figure 2(a) shows the sum 
of the photoelectric detection event probabilities at the 
two detectors P = p\ +P2 as a function of the single-atom 
Rabi angle 9, varied by changing the intensity of the fl± 
field. The data are fit with the sinusoidal oscillation of 
Eq. (1) modified by two Gaussians, as described in the 
Methods section. The choice of the fit function is moti- 
vated by a physical picture in which the visibility of the 
oscillation is smeared by fluctuations of the atom num- 
ber and the intensities of laser fields fii and fi 2 , while the 
leading contribution to the overall decay of the retrieved 
signal, for the case when y^Ve^o < AEb, is due to an in- 
homogeneous distribution of light shifts for atoms in state 
\R), ~ N c £Iq/AE b which couple the state \R) to other 
singly-excited states. The effective number of atoms N e 
is defined as N e = ^j=i /^ 2 (0)- For our experi- 



mental geometry flf = fi 2 (0) exp(— 2x 2 /w 2 — 2y 2 /w 2 ), 
and the atom density p = po exp(— 2z 2 /w1). Therefore 
7V e = {-K/2f/ 2 w x w y w z p . 

To explore the collective character of the observed Rabi 
flopping, we measure P as a function of 9 while varying 
the peak density of the sample po, Fig. 2(b-d). Figure 
2(e) shows the normalized frequency of the Rabi oscilla- 
tion ^/Wl extracted from the data in Figure 2(a-d) as a 
function of number of atoms in the ensemble N a . The 
latter is calculated using peak density po measured by 
the hyperfine state-selective fluorescence imaging of the 
atomic sample with magneto-optical trap cooling beams 
used without a repumping field to exclude contribution 
of \5si/2,F = 1) atoms. The absence of additional peaks 
in Fig. 1(c) supports a near-unity value for the fraction 
of atoms / in the m = Zeeman sub-level. Ideally, we 
expect the effective atom number N e extracted from the 
Rabi oscillation period to be equal to the atom number 
N a determined by the fluorescence imaging of the sam- 
ple, so that parameter C in the fit in Fig. 2(e) would 
equal unity, whereas we extract C = 0.74. In addition 
to the factor y/J, likely cause for C < 1 are alignment 
imperfections, uncertainties in the determined waists of 
the two-photon excitation laser beams, and uncertainties 
in the fluorescence measurements of po. 

We further confirm that the dynamics seen in Figure 
2 correspond to the oscillation of Eq. (1) by measure- 
ments of the second-order intensity correlation function 
at zero delay (0) as a function of 9, shown in Figure 3. 
Measured values of (0) well below unity, together with 
substantial visibility of the oscillations, indicate that only 
one Rydberg excitation is present in the entire ensemble 
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FIG. 5: Probability of photoelectric detection P as a func- 
tion of the single-atom Rabi angle 9. Level |102si/ 2 ) is excited 
for r = 0.2 /is in a, and level ]102d 3 / 2 ) is excited for r = 1 /is in 
b. The solid curves are fits as in Figure 2(a-d), the fit parame- 
ters (A,a,b,N e ) are (4.56,5.27,3.86,340) in a. The excitation 
spectrum for |102d 3 / 2 } shows a complex structure due to an in- 
terplay of the bias magnetic field and an ambient electric field. 
For the data in b, the laser is tuned to the strongest spectral 
component, with the scale 9' determined by using the value of 
N e = 492 from the preceding measurements with the |102si/2) 
level, with a fit providing the value of peak single-atom Rabi fre- 
quency flo and the fit parameters (A, a, b) are (2.58, 10.7, 3.49). 
The vertical error bars represent ± one standard deviation (vM) 
for M photoelectric counting events. The horizontal error bars 
in b reflect the uncertainty in determination of the x-axis scale 
9'. 

of several hundred atoms. The measured values of g( 2 > (0) 
for \/Wl6 > 57r in Figure 3(b) suggest that the contribu- 
tion of populated doubly-excited states to extracted val- 
ues of a is substantial. On the other hand, the sizable os- 
cillation visibilities in Figure 2 indicate small population 
of the doubly-excited states for \/Wl9 « 7r, as the lat- 
ter are expected to have a suppressed retrieval efficiency. 
Combining all the data points for \/WlQ ~ it in Fig. 3(b) 
we obtain g^(0) = 0.006(6), which to our knowledge is 
the lowest value for this quantity for any previously re- 
ported light source. It is consistent with a lower bound 
of gl g (0) = 0.012(2) due to background counts, of which 
about half are due to detector dark counts. 

The importance of the condition AEb 3> Slj for obser- 
vation of many-body Rabi oscillations is checked by re- 
ducing AEb in measurements with n = 90 and n = 81, 
Fig. 4. Figure 5(a) shows data with increased 5u by 
using shorter r = 0.2 /is excitation. The oscillation visi- 
bility is clearly lower both for smaller AEb, Figure 4, and 
larger 5u), Figure 5(a). Figure 5(b) shows an oscillation 
with a similarly reduced visibility in measurements with 
the |100g? 3 / 2 ) level with r = 1 /jjs excitation, which may 
be attributed to the blockade breakdown due to a strong 
angular dependence of the atomic interaction strengths 



for | nd) -levels [l3j . It should be also noted that for a 
Gaussian distribution of atom-field couplings, describing 
a laser beam exciting the extended gas in our experi- 
ment, the single-atom Rabi oscillations are almost com- 
pletely washed out (23|, which makes the observation of 
many-atom oscillations under these conditions even more 
remarkable. 

We have demonstrated coherent many-body Rabi oscil- 
lations in an ensemble of several hundred cold rubidium 
atoms. The oscillations provide compelling evidence for 
the achievement of a collective Rydberg excitation block- 
ade by a single excited atom. Our results pave the way 
for the realization of quantum simulators which employ 
qubits encoded in atomic ensembles, with fast two-qubit 
quantum gates mediated by strong Rydberg interactions 

In. 

METHODS 

Atom preparation and laser excitation. A magneto-optical trap 
of 87 Rb is loaded from background vapor for 70 ms. During the fol- 
lowing 25 ms, the detuning of cooling light is increased, the rcpumpcr 
intensity is decreased, and the optical lattice is turned on. The lat- 
tice is composed of a single 782 nm retro-reflcctcd linearly polarized 
Gaussian beam. Untrappcd atoms arc allowed to fall away from the 
experimental region during the next 15 ms period, and a Bq — 4.3 
G bias magnetic field is turned on. The trapped atoms arc optically 
pumped to the |5s 1 / 2 i^ 7 ' — 2,mjj- — 0) state, the optical lattice is 
switched off by an acousto-optical modulator (AOM), and a 3 fis long 
sequence of two-photon Rabi driving and retrieval is repeated for 50 
/is, with a 1 fis optical pumping period included every five cycles. The 
overall repetition rate of the experiment is ~ 8 Hz. For the data in 
Figure 2 the peak density po was controlled by varying the time pe- 
riod between lattice loading and the two-photon excitation sequence 
between 15 and 90 ms. Both Q± (at 795 nm) and H2 (at 475 nm) fields 
arc linearly polarized along the same axis. The 795 nm field is produced 
by an extended cavity diode laser (ECDL). Light at 475 nm is gener- 
ated by frequency-doubling the output of a tapered amplifier driven 
by a 950 nm ECDL laser. Both lasers are frequency-locked to a ther- 
mally stabilized ultra-low expansion glass cavity and have linewidths 
< 100 kHz. The transition is located by scanning the laser frequency 
across a resonance and measuring the photoelectric detection probabil- 
ity for the retrieved field. The 795 nm and 475 nm excitation fields arc 
tuned to the two-photon resonance between the ground-level compo- 
nent |5si/2j F = 2, mp — 0) and a Zccman component of the Rydberg 
level |ns 1 / 2 i m j — — 1/2). The single-photon Rabi frequency on the 
blue transition is Q2 — — e £{5pi/2 F — 1» irij? — 0|r|nsi /2, m i) > where 
S is electric field amplitude. The radial matrix clement is reduced 
using the Wigncr-Eckart theorem. The angular part is calculated fol- 
lowing Ref. |25H , while the reduced matrix element is approximated by 
(r) = 0.14 x (50/n) 3/2 a fll . 

Since Q2 and fields arc propagating in the same spatial mode, the 
retrieved field is phase matched into the mode of the Oi field and cou- 
pled into a single mode 50/50 fiber beamsplitter followed by a pair of 
single-photon detectors Di and D2. A gating AOM at the fiber beam- 
splitter input port is employed to avoid damaging the single photon 
detectors by the f2i field. 

Photoelectric detection and data analysis. For every experi- 
mental trial, photoelectric events on detectors Di and D2 arc recorded 
within a time interval determined by the length of the retrieved pulse. 
Photoelectric detection probabilities for both detectors are calculated 
as pi,2 — iVi^/iVfl, where N± } 2 are numbers of recorded events and Nq 
is the number of received triggers. The photoelectric detection prob- 
ability for double coincidences is calculated as N12/N, where 7Vi2 is 
a total number of simultaneous clicks on both detectors for a given 
experimental trial. The second order intensity correlation function at 
zero delay is given by g 2 (0) — P12/ (piPz)- Transmission through the 
glass vacuum chamber is 0.92, the gating AOM diffraction efficiency 
is 0.7, the fiber coupling efficiency is 0.73, and the quantum efficiency 
of the single-photon counters is 0.55. Normalizing the probability of 
photoelectric detection event per trial P{9 tt / V N e ) 0± 0.04 by the 
product of these efficiencies we arrive at the single-photon generation 
efficiency into the Gaussian mode of our single-mode fiber e — 0.15. 

Decoherence model. We employ the following Hamilto- 
nian to describe our system: H — h{uj g a^ L 9 + UJ r cr£ r ) + 

i ft(n M e~ iw *'*££ fl + h.c.) + J2n>v hA ^ a 7 ® K r - Thc atomic op- 
erators for the atom \x are defined as <r^ b — |a)^{6|, where a, b £ [g, r] 
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with \g) fj, being the atomic ground state and \r) ^ being the addressed 
Rydberg level. The two-photon excitation is modeled using the effective 
Rabi frequency O — Hi 0,2 / (2 A) . The interaction between Rydberg lev- 
els is described with a single-channel model. For A ML , O^, V(/i, f ), 
the excitation blockade is operational. Adiabatic elimination of double 
and higher-order excitations from the equations of motion results in 
an effective Hamiltonian for the singly-excited part of the spectrum: 

Heft = E 3 ^MU\+T, i>j nc i:l (\i)ti\ + \j)(i\) + h £j &%(U>(G| + 

\G)(j\). Here Aj = - J2i^j n i/( 4A * 3 ). Cij = -O i f2 j /(4A ij ), where 
\j) is the many-body state with the j-th atom in the Rydberg level. 
The first two terms of the effective Hamiltonian are due to the light 
shifts induced by the (detuned) doubly-excited states onto the single 
excitations. 

When the interaction-induced inhomogencous light shifts are 
omitted, the Hamiltonian results in an ideal Rabi oscillation be- 
tween the ground state \G) and the single spin wave \R) — 



(1/y O 2 ) Y2j If a t ti me t — the system is in state |G), the 

state at future times is given by |^(r.)) — cos (Or /2) | G) —i sin (Or /2) | R) . 
When the light shift terms are included, the state \R) is coupled to 
a broad distribution of singly-excited states and therefore leaks into 
this quasi-continuum, leading to P ~ | (i?|'0(r.)) | 2 decaying with a rate 
~ N c Qq/AEb- The doubly-excited states are expected to be populated 
at a rate ~ 7V e 2 . Trial-to-trial fluctuations AO and AN e in N e and 
Oo, respectively, lead to a decay of the oscillation visibility. The proba- 
bility of photoelectric detection per trial P as a function of 9 in Figures 

2-5 is therefore fit by a function ^Ae~~ a() (l — e~ b ® cos(V N e 8)), where 
dimcnsionless fit parameters a ~ iV e and b ~ (AO/Oo) 2 -t- (AN e /2N e ) 2 
describe the roles of the light shifts and population of doubly-excited 
states, and atom number and intensity fluctuations, respectively, while 
an amplitude A represents the overall measured retrieval and detection 
efficiency. 
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